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ABSTRACT 


The  feasibility  of  utilizing  an  atmospheric  velocity 
gradient  for  obtaining  continuous  non-powered  flight  was 
analyzed  mathematically  and  experimentally.  It  was  dis- 
covered that  there  seems  to  be  no  theoretical  bar  to  using 
wind  shear  for  this  purpose.  The  velocity  gradients  in 
the  upper  atmosphere  were  found  to  be  too  weak  for  main- 
taining continuous  non-powered  flight,  but  in  the  vicinity 
of  the  earth's  boundary  layer  wind  shear  strengths  nec- 
essary for  this  purpose  were  observed. 

A theoretical  flight  pattern  was  developed,  numerically 
analyzed  and  flight  tested  using  a powered  aircraft  to  sim- 
ulate high  performance  sailplanes.  Though  an  increase  in 
aircraft  performance  was  realized  a zero  net  loss  of  altitude 
was  not  achieved.  It  is  concluded,  however,  that  further 
research  is  warranted,  especially  concerning  the  related 
topics  of  jet  stream  shears,  vertical  updraft  considerations 
and  more  efficient  flight  schedules. 
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INTRODUCTION 

Numerous  articles  have  been  written  dealing  with  the 

meteorological  phenomenon  known  as  wind  shear.  Wind  shear 

is  a velocity  gradient  with  increasing  altitude i a wind  pro- 

file  which  must  exist  in  order  to  account  for  the  difference 

between  relatively  large  upper  atmosphere  winds  and  the  air 

at  the  earth's  surface  which  is  essentially  at  rest. 

Such  a wind  differential  has  become  an  important  design 

consideration  of  modern  skyscrapers.  Structural  allowances 
* 

must  be  made  for  the  increased  wind  loading  experienced  by 

. 

tall  buildings.1  Recently  studies  concerning  the  effects  of 
velocity  gradients  upon  landing  aircraft  have  been  undertaken. 

The  National  Aeronautics  and  Space  Administration  sponsored 
an  investigation  into  such  effects.2  Analytic  studies  util- 
lzing  a digital  computer  simulation  show  that  wind  shear  pro- 
duces a significant  deviation  in  the  touchdown  point  of  a land- 
ing aircraft  as  compared  to  the  same  aircraft  landing  under 
exactly  the  same  conditions  except  for  the  absence  of  shear, 3 
Additionally,  the  National  Transportation  Safety  Board  acknowl- 
edged wind  shear  as  a' primary  contributor  to  the  crash  of  an 
Iberia  airliner  on  December  1?,  1973.4  Therefore  the  exist- 
ence of  velocity  gradients  has  emerged  in  recent  years  as  an 

• I 
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important  engineering  consideration. 

Types  Of  Wind  Shear 

Evidence  suggests  that  there  exist  two  characteristic 
types  of  wind  shear.  The  first  is  found  near  the  earth's 
stirface  and  closely  resembles  the  shape  of  a typical  viscous 
flow  boundary  layer.  Beyond  the  well  defined  boundary  layer, 
however,  is  a more  subtle,  gradually  increasing  velocity 
gradient.  This  gradient  has  been  referred  to  as  upper  level 
shear. 

Boundary  Layer  Shear 

A great  deal  of  research  has  been  done  on  the  subject 
of  the  earth's  boundary  layer  shear.  This  is  probably  due 
to  the  previously  mentioned  importance  of  boundarj  lajor  sheai 
in  recent  engineering  considerations  and  its  relative  ease 
of  measurement  due  to  its  close  proximity  to  the  ground. 

Boundary  layer  shear  is  a relatively  common  occurrence 
because  of  its  dependanee  upon  surface  friction  effects  rather 
than  prevailing  weather  conditions.  The  single  major  exception 
to  this  generalization  is  pointed  out  by  Mr.  Davenport  in  his 
initial  study  of  boundary  layer  shear.  He  found  that  in  severe 
thunderstorms  or  squalls  the  boundary  layer  breaks  down.  Under 
such  violent  frontal  conditions  the  atmosphere  directly  above 
the  earth  is  highly  turbulent  and  unstable,  the  result  being 
that  a distinct,  gradient  style  increase  in  velocity  with  al- 
titude is  greatly  reduced. 5 
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Upper  level  Shear 

Analytic  research  concerning  upper  level  shear  has  been 
much  less  extensive  than  for  boundary  layer  shear.  High  al- 
titude shear  requires  much  more  sophisticated  measuring  sys- 
tems and  is  of  less  direct  engineering  importance  due  to  its 
relative  weakness  and  its  distance  from  the  earth's  surface. 

Actual  meteorological  data  confirms  the  existence  of  upper 
level  wind  shear.  Figure  2 displays  a velocity  gradient  which 
occurred  over  Omaha  in  January,  1974.  The  strength  of  the 
gradient  was  ,0065  ft/sec^ft,?  Two  important  characteristics 
may  be  noted  about  the  non-boundary  layer  shear  described  in 
Figure  2,  First,  as  mentioned  previously,  the  gradient  is 
quite  weak,  increasing  only  .65  ft/sec  in  every  100  feet. 

This  fact  alone  suggests  that  detection  and  accurate  measure- 
ment of  upper  level  shear  would  demand  painstaking  care.  Sec- 
ondly, the  velocity  gradient  above  the  boundary  layer  may 
reasonably  be  represented  by  a straight  line  function. 

In  order  to  achieve  maximum  efficiency  from  flying  time 
for  this  report  it  was  necessary  to  determine  exactly  what 
type  of  weather  conditions  vrould  prove  most  advantageous  for 
detecting  and  measuring  the  velocity  gradient.  In  particular, 
a forecast  outline  for  wind  shear  was  needed  for  the  Virginia 
and  Maryland  area  where  data  flights  were  to  be  made, 

Mr,  W,  Bonner  discussed  the  characteristics  of  low  level 
wind  maxima  in  an  article  written  for  Monthly  Weather  Review. 


n 
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Though  the  thrust  of  the  article  was  aimed  at  commenting  on 
characteristics  of  the  boundary  layer  particular  to  the  mid- 
•west,  the  author  made  valuable  comments  concerning  what  might 
be  expected  of  a velocity  gradient  on  the  Eastern  seaboard. 

Also,  it  is  reasonable  to  assume  that  where  and  when  a pro- 
nounced boundary  layer  shear  exists,  a corresponding  upper 
level  gradient  may  appear.  Therefore  Mr.  Bonner's  analysis 
of  weather  conditions  creating  low  level  shear  may  also  be 
applied  to  predicting  upper  level  shear.  His  study  was  based 
on  data  taken  from  forty-seven  weather  stations  (including 
several  in  the  Maryland,  Virginia  and  Delaware  areas)  during 
a two  year  time  span  at  a rate  of  two  observations  per  day. 

In  general,  the  velocity  gradients  were  more  frequently 
observed  in  the  morning  hours.  Though  in  the  midwest  there 
was  a definite  increase  in  low  level  shear  during  the  summer 
months,  no  such  precise  generalization  co\ild  be  made  con- 
cerning the  East,  Actually,  the  velocity  gradient  seemed  to 
be  more  abundant  in  the  Annapolis  area  between  the  months  of 
October  and  March,  Tn  addition,  such  a velocity  gradient  v*s 

O 

usually  from  the  north  or  northwest. 

An  independant  U.  S,  Army  Missile  Command  report  supports 
Mr,  Bonner's  observations.  It  shows  an  increase  in  the  strength 
of  the  shear  in  the  vicinity  of  Annapolis's  longitude  during  the 
fall  and  winter  month  s,9  Figure  3 Is  a graphical  display  of 

I 

some  of  the  results  of  this  report. 


APR  MAY  JUN  XL  AUG  SEP  OCT  NOV  DEC 
TIME 


Figure  3*  Time  section  of  wind  shear 
(m/see»km)  near  longitude 
80°W,  Annapolis  longitude 
is  approximately  76°30'W. 

Based  upon  these  two  reports  the  following  guidelines 
were  devised  to  govern  the  choice  of  days  for  data  collecting 
flights. 

If 

Flights  should  me  made? 

a. )  during  the  late  fall,  winter  and  possibly  early 
spring  months. 

b. )  over  the  ocean,  similar  large  body  of  water  or 
above  about  one  thousand  feet  to  avoid  ground  induced 
turbulence  and  friction  effects. 

c. )  during  the  morning  hours  to  avoid  afternoon  con- 
vection currents. 


d.)  on  a day  with  strong  northerly  surface  winds  fol- 
lowing the  passage  of  a low  pressure  cold  front. 
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ANALYTIC  DISCUSSION 


The  effect  of  wind  shear  upon  an  aircraft  is  that 
of  either  acceleration  or  deceleration,  depending  upon 
whether  the  speed  in  inertial  space  (ground  speed)  is 
Increasing  or  decreasing.  Since  the  shear  is  a vel- 
ocity differential  with  altitude  the  aircraft  must  not 
be  moving  perpendicular  to  the  wind  direction.  In 
crosswind  flight  no  longitudinal  acceleration  or  de- 
celeration due  to  the  wind  shear  would  occur. 

Figure  4 shows  the  aircraft  center  of  gravity  de- 
scending in  a wind  shear  parallel  to  the  wind  direction. 
The  positive  directions  have  been  indicated.  The  flight 
angle  If  has  been  chosen  positive  clockwise  from  the  x 
axis.  Summing  the  forces  in  the  x and  y directions i 

£ Fx=Max=Lsin)f-DcosX  ( \ ) 

£Fy=0=Lcosl(-W+  Dsintf  (Dsin  t ~ 0)  (2) 

The  wind  direction  has  been  assumed  to  be  unidi- 
rectional, acting  only  parallel  to  the  horizontal  axis, 
and  the  gradient  shape  represented  by  a straight  line 
function.  Since  the  wind  direction  is  assumed  to  act 
only  parallel  to  the  x axis,  only  inertial  acceleration  in 
that  direction  may  exist. 
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Figure  4.  Aircraft  center  of  gravity  descending  with 
linear  shear. 

Small  angle  approximations  may  be  applied  because  the 
descent  angle  for  best  glide  flight  is  in  general  quite 
small.  This  is  particularly  true  for  high  performance  sail- 
planes and  even  holds  true  for  most  general  purpose  pow- 
ered aircraft. 

Therefore;  costf=l 

sin?=Y 

tan*=l( 

where  <S  is  measured  in  radians. 

Equations  (1)  and  (?)  may  now  be  simplified; 


Max=lJ-D  (3) 

L=W  (4) 

Combining  equations  (3)  and  (4)  yields 
Max“WJf-D 

' whore  ^qSfC^in+Cy/VtrAt  e)  (6) 
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The  aircraft's  groundspeed  is  the  sum  of  its  true 
airspeed  and  the  wind  speed; 

Vg=V+Vw  (9) 

Here  it  may  be  noted  that  again  the  small  angle 
approximation  cosX  has  been  applied  so  that  V=VcosJf. 

If  the  velocity  is  held  constant  then  the  acceleration 
becomes 


* ax~dVw 
dt 


(10) 


The  change  in  wind  speed,  however,  may  be  considered 
the  product  of  the  rate  of  descent  and  the  velocity  gra- 
dient; 


dVws-dVw , dh 
dt  dh  dt 


(ID 


This  product  is  negative  (dh/dtcO)  indicating  that 
the  acceleration  is  in  fact  a deceleration;  the  airplane 
is  losing  ground  speed  while  maintaining  a constant  air- 
speed. 

Equation  (11)  may  be  further  rewritten  realising  that 
the  rate  of  descent  dh/dt  is  also  the  downward  velocity 
Vy=VsinV, 








dVw  = -dVw  Vy  = -dVw  V* 


Combining  equations  (5),  (.6)  and  (12)  yields 

-WVjfdVw  = W#  - qSC^in  - qC^S/ft/Re 
B dh 

but  CL2  = (L/qS)2  = (W/qS)2 

so  equation  (13)  becomes 

-WlfdVw  = W8  - qSCDmin  - W^qSrr/R  e 
B dh 

Dividing  through  by  w/g  and  rearranging  terms; 

-VVdVw  - g*  = -g/qSCpHin  + w ^ 

1111  1 W qStrfi  <J 

Solving  for  the  angle  of  descent; 


(g  + Vd Vw) 
dh 


AscDmin  + 

V w 


W 

qSnJR  e 


With  no  velocity  gradient  present  (dVw/dh=0  or  cross- 
wind  flight); 

Y = + w = D 

W qSH/R  e L 

Equations  (17)  and  (18)  show  that  a constant  velocity 
descent  with  a wind  shear  differs  from  a descent  without 
a shear  (or  crosswind  to  a shear)  by  a factor  F; 

(g  -f  VdVw)  ^ 

dh 

Taking  as  an  example  the  shear  of  Figure  2 of  the 
Introduction,  with  a velocity  of  100  mph  = 146.6?  ft/sec, 

F = .97.  Thus  a flight  angle  is  reduced  three  percent  in 
a wind  shear  of  .006$  ft/sec/ft  at  100  mph.  Using  a stronger 
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shear  of  .02  results  In  F = .92,  an  eight  percent  reduction 
of  the  flight  angle.  A reduction  of  flight  angle  (in- 
• crease  in  F)  nay  also  apparently  result  from  an  increase 
in  true  airspeed  V,  However,  it  must  be  remembered  that 
this  analysis  has  at  its  foundation  the  small  angle  ap- 
proximations which  may  be  violated  at  higher  airspeeds  in 
the  descent. 

Equation  (l6)  provides  a theoretical  method  for  meas- 
uring a straight  line  wind  shear j 

<*» 

where  R/D  is  the  rate  of  descent  at  constant  velocity  V. 

R/D  = h2  - hj 
*2  - tl 

Since  q depends  upon  temperature  and  the  aircraft 
weight  will  vary  due  to  the  loss  of  fuel  (powered  air- 
craft assumed)  an  average  of  those  two  quantities  must  be 
determined  for  a given  descent.  An  overly  long  descent 
might  cause  errors  due  to  inaccurate  temperature  or  weight 


averages  being  applied  to  equation  (20). 

Though  the  foregoing  analysis  demonstrates  the  theo- 
retical effect  of  wind  shear  upon  an  aircraft,  it  is  of 
little  practical  value.  Glide  angle  reduction  in  a con- 
stant velocity  descent  is  quite  small.  To  reach  a sig- 
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shear  would  bo  required.  Shears  of  that  order  of  mag- 
nitude have  not  been  observed  in  nature.  Also,  though 
a similar  interpretation  of  a constant  velocity  climb, 
into  the  wind  yields  a result  equivalent  to  the  descent, 
it  is  obvious  that  a constant  velocity  climb  is  not 
practical  for  any  non- powered  aircraft.  Another  problem 
with  this  analysis  is  that  it  does  not  allow  for  a con- 
venient interpretation  of  the  transition  from  the  descent 
to  the  climb.  Some  type  of  realistic  climb  and  descent 
schedule  must  be  devised  which  can  be  followed  by  the  pilot 
without  excessive  difficulty. 

Consider  a purely  circular  flight  path.  Figure  5 
shows  t.ho  aircraft's  movement  as  projected  onto  the  ground. 
Flight  path  heading  is  measured  clockwise  from  the  wind 
direction  to  the  aircraft  velocity  vector. 
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The  aircraft’s  groundspeed  is  the  sum  of  its  true  air- 
speed and  the  component  of  wind  acting  along  the  flight 
path  | 

Vg  = V + Vwcos V (21 

The  inertial  acceleration  experienced  by  the  air- 
craft due  to  climbing  and  descending  through  the  shear 
is  again  the  change  in  groundspeed. 

ax  = dVg  = dV  - VwsinVd'f  + cos^dVw  (2J 

dt  3t  dt  dt 

where  the  x direction  is  now  taken  to  correspond  to  the 

direction  of  the  flight  path. 

For  this  discussion  the  velocity  of  the  aircraft  is 
limited  to  a sinusoidal  variation  so  over  one  period  of 
the  flight  cycle  its  derivative  may  bo  assumed  zero. 

Altitude,  true  airspeed  and  climb  angle  may  be 
written  as  functions  of  the  flight  angle.  The  climb 
angle  will  be  zero  at  T = 270* and  90*  . It  will  reach 
a maximum  positive  value  at  V =180*  and  a minumum  at 
V = 0° . Velocity  and  altitude  will  oscillate  about 
median  values  V and  h respectively  in  a sinusoidal  manner. 
Airspeed  reaches  a maximum  at  V =9$ while  the  aircraft 
is  at  the  highest  point  in  the  cycle  at  V = 2 70° . These 
quantities  may  be  converted  to  the  time  domain  by  assum- 
ing a constant  turning  rateil. 

M*  =nt 


i 


(23a) 


dY  =iidt 


(23b) 


The  equations  for  altitude,  velocity  and  glide  angle  are; 

h = h - ^AhsinY  (2 

V = V + | AVsinY  (2 

^ ~ “)fmaxcos  (2 

where 

V = ^Vmax  + vmin)/2 

AV  = <Vmax  - vmin) 

Ah  = (hmax  - hpj^y,) 

Returning  now  to  equation  (22),  substituting  in 
equation  (11)  for  the  rate  of  change  of  the  wind  vel- 
ocity and  multiplying  by  the  aircraft  mass  produces  the 
inertial  fo^ce  due  to  the  wind  shear  acting  on  the  air- 
craft s 

F «=  W /cosYdVw  dh  - VwsinYdyN  (2 

g l dh  dt  dt J 

The  term  F in  this  equation  should  not  be  confused 
vith  the  factor  determined  in  equation  (19).  When  this 
force  is  less  than  zero,  the  airplane  is  decelerating 
in  inertial  space.  This  occurs  during  the  climb  vihen 
the  flight  heading  is  between  90° and  2?0° , During  this 
interval  the  aircraft  is  extracting  energy  from  the  wind. 
Similarly  in  the  descent  the  aircraft  is  accelerating 
relative  to  the  ground  and  it  must  give  up  energy  by 


• either  losing  altitude  or  by  providing  thrust  from  an 


engine. 
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The  power  available  to  the  aircraft  due  to  the  wind 
shear  is; 

Pa  = -FV  a -VW  /cosYdVw  dh  - VvsinY_fL\  (2 

g V 3K”  dt  J 

Negative  power  represents  power  supplied  by  the  wind 
shear.  Equation  (23b)  has  been  substituted  into  equation 
(28)  to  replace  the  time  rate  of  change  of  the  aircraft 
flight  heading. 

The  climb  or  descent  rate  (dh/dt)  must  now  be  ana- 
lyzed in  terms  of  the  flight  heading.  Assuming  the  ex- 
change of  energy  during  a climb  or  descent  may  be  repre- 
sented by  a loss-gain  relationship  between  kinetic  and 
potential  energy; 

4PE  =«KE 

mgAh  = ImAV2 

Ah  = AV2 /2g  = (y&ax  - v£in)/2g  (2 

also i the  climb  rate  will  vary  throughout  the  cycle, 


dh/dt  = -(dh/dt)M  cosy 

The  subscript  M Indicates  maximum  climb  rate.  Using 
equation  (23b)  and  integrating  over  the  climb  half  of  the 
period; 

(270°  (270* 

Ah  a \ dh  = -(dh/dt)M  \ cosfdV 
' 90®  —rr^JW 

Ah  = 2(dh/dt)M  f7 

__ 


'■'qua t. ions  (?Ql  and  (30  may  be  combined  t.o  discover 
the  maximum  rate  of  climb; 

(dh/dt)M  =_aaV?/4g  ^2 

Now  with  the  exception  of  the  wind  speed  the  power 
available  may  be  written  in  terms  of  the  flight  heading 
by  substituting  equations  (25),  (30)  and  (3?)  into  equa- 
tion (?R) ; 

Pa  = W£L(V  + 2 flVsinf ) (cos^Y  dVw  AV^  + VwsinV)  (33) 

f dh  4g 

The  wind  speed  and  its  derivative  are  also  functions 
of  the  flight  heading  and  vary  in  a sinusoidal  manner 
as  the  aircraft  traverses  the  flight  path. 

Both  boundary  layer  shear  and  upper  level  shear  may 
be  represented  by  the  same  power  equation; 

Vw  = kh1^  = k(h  - ^Ahsin  )X  (34a) 

dVw/dh  = k h1/a"t  = k(h  + +4hsin  )' ^a_1  (34b) 

a a 

For  boundary  layer  shear  the  Davenport  parameters 
discussed  previously  may  be  applied  while  if  a linear 
shear  is  experienced  a = 1 and  k is  the  slope  of  the  shear. 

Now  an  equation  for  power  available  may  be  written 
completely  in  terms  of  the  flight  heading. 

Pa=^V+^AVsinY)(cos?m^  + k(h-iAhsinV)sinY)  (35) 

g 4g 


A linear  shear  of  slope  k has  been  assumed  for  equa- 
tion (35).  Consequently  h may  need  ad.iustment  so  that  the 
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shear  graph  Intersects  the  horizontal  axis  at  the  origin. 
(Note  the  extrapolation  of  Figure  2 to  the  velocity  axis.) 
In  equation  (35)  k = dVw/dh. 

Equations  (31)  and  (32)  may  now  be  combined  to  re- 
place iih  in  the  power  equation!  v 

Pa  * WTl  (V+lAVsinV)/ kAV2cos2y  +khsinY  -kA^sin2^  ) 

g \W  ) 

but  (eos2y-  sin2^  ) = eos2V  . 

Therefore  the  power  equation  finally  reduces  to 

Pa  » Wil(V  + iAVsinV )(kAV2eos2y  + khsinY)  (' 

g 

At  this  point  it  is  desireable  to  determine  the 
energy  provided  to  the  aircraft  due  to  the  wind  shear  in 
one  cycle  and  compare  with  the  energy  required  (drag 
effects)  during  the  course  of  the  same  cycle.  The  two  re- 
sulting energy  expressions  may  be  equated  to  determine  the 
parameters  necessary  to  prevent  a loss  of  altitude  in  the 


cycle. 


The  energy  available  to  the  aircraft  due  to  the  wind 


shear  over  one  cycle  is! 

Ea  a j Pa  dt  = 1 Jpa  dV 


Integrating  equation  (3^)  produces  an  expression  for 
the  energy  available! 


WAV  IT  f V2  k + Vw  j 
2*  \ g ) 


i 


The  energy  required  to  complete  one  cycle  is  the 
integration  of  the  power  required  expression  of  the  air- 


?0 


] 

. 


i 
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craft. 

Flight  testing  of  tha  circular  flight  path  proved  un- 
satisfactory and  impractical  even  under  simulated  high 
glide  ratio  conditions  (see  Appendix  1).  The  following 
observed  conditions  contributed  to  the  inadequacy  of  the 
circular  approach  to  the  dynamic  soaring  problem. 

The  flight  schedule  demanded  far  too  much  effort 
from  the  pilot  in  order  to  conform  exactly  to  the  the- 
oretical model.  Coordination  of  flight  heading,  velocity 
and  altitude  proved  to  be  a formidable  task.  It  was  vir- 
tually impossible  to  obtain  the  exact  sinusoidal  var- 
iations of  the  flight  parameters  as  required  by  theory. 

The  result  was  that  a truly  constant  turning  rate  was  not 
achieved}  bank  angle  varied  from  a spectacular  nearly 
vertical  wing-over  at.  the  high  altitude  point  and  a med- 
ium hanked  turn  at.  the  low  altitude  end.  Also,  as  may 
be  ascertained  from  Figure  *5,  the  aircraft  was  not  taking 
maximum  advantage  of  the  wind  shear  effect  through  most 
of  the  cycle.  The  component  of  wind  acting  along  the  lon- 
gitudinal axis  of  the  aircraft  varies  sinusoidally.  Thus 
the  flight  path  obviously  reduces  the  effect  of  the  wind 
shear  upon  the  flight  angle  in  the  ascent,  or  descent  phases 


I 


of  the  flight  schedule 


From  these  unsatisfactory  characteristics  arose  the 
desire  to  modify  the  circular  flight  path  by  making  the 
ascent  and  descent  legs  entirely  parallel  to  the  wind  dir- 
ection. Figure  6 depicts  such  a flight  schedule  and  shall 
be  hereafter  referred  to  as  the  "racetrack"  flight  path. 

The  racetrack  flight  schedule  may  be  broken  into 
three  portions  for  examination.  The  first  is  the  com- 
bination of  the  ascent  and  descent  legs,  the  second  is 
the  high  speed  turn  at  the  low  altitude  end  and  the  third 
is  the  low  speed  turn  at  the  high  altitude  end. 

In  order  to  maximise  efficiency  of  the  aircraft's 
exchange  of  potential  and  kinetic  energy  in  the  climb 
and  descent,  the  phugoid  oscillation  mode  of  the  aircraft 
will  be  utilized.  The  phugoid  oscillation  is  a stick- 
fixed,  long  period,  lightly  damped  natural  longitudinal 
oscillitory  mode  in  which  the  aircraft  alternately  ex- 
changes airspeed  for  altitude  at  a natural  frequency. 

A detailed  mathematical  evaluation  and  discussion  of  the 
phugoid  oscillation  is  contained  in  Etkin's  Dynamics  of 
Flight  or  other  similar  stability  and  control  texts. 

The  velocity  of  the  aircraft  in  the  ascent  and  descent 
legs  is  described  by  the  relation 

V = V + AV/2  coswt  (39 


Figure  6.  Racetrack  1 light  schedule. 
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V and  AV  possess  the  sane  definitions  as  in  the  circular 
case,  w is  the  phugoid  angular  frequency  and  f is  the 
phugoid  frequency  in  cycles  per  second. 

The  sign  of  the  cosine  term  in  equation  (39)  depends 
upon  where  aero  time  is  established.  If  zero  time  is  con- 
sidered to  be  at  the  high  altitude  low  speed  end  the  sign 
is  negative.  Conversely  if  t*0  at  the  high  speed  end  the 
sign  is  positive. 

At  each  end  of  the  phugoid  oscillation  a turn  at  con- 
stant altitude  and  velocity  is  assumed.  Each  turn  must 
be  tailored  for  minimum  energy  expenditure  while  they  must 
also  necessarily  be  limited  by  phugoid  maximum  and  min- 
imum velocities,  structural  limitations  (as  determined  by 
the  manufacturer's  recommended  maximum  load  factor),  and 
pilot  physiological  response  to  turn  loading. 

The  advantages  of  this  flight  path  as  opposed  to  the 
circular  path  are  readily  discerned.  The  ascent  and  de- 
scent legs  are  performed  entirely  parallel  to  the  wind 
direction  and  therefore  make  maximum  use  of  the  longitudinal 
wind  shear  effect.  Also,  this  schedule  makes  it  easier  for 
the  pilot  to  conform  to  theory. 

Piloting  procedure  is  merely  toj 

1.  Create  the  phugoid  oscillation  with  the  desired 
velocity  parameters. 


2,  Provide  a "tap"  to  the  stick  in  the  descent  to 
counteract  the  natural  phugoid  damping  and  maintain  vel- 
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ocity  parameters, 

3.  At  each  extreme  altitude  point  in  the  flight 
schedule  complete  a constant  velocity  turn  of  1B0  degrees 
and  release  the  stick  so  that  the  aircraft  may  resume  its 
oscillation. 

The  power  required  for  the  aircraft  to  complete  one 
cycle  of  the  racetrack  flight  schedule  is  the  sum  of  the 
power  required  in  both  turns  and  the  power  required  for 
one  complete  phugoid  oscillation. 

The  power  required  for  a level  turn  isj 

rr---rv,*(CDwln  + nV  ) (41) 

pW4 

For  convenience  the  constants  ± S.  »«d  /TT  « ha'*c 
been  replaced  with  P and  A respectively. 

The  minimum  drag  coefficient  is  found  by  taking  the 
derivative  of  the  power  with  respect  to  velocity  and  set- 
ting the  resulting  expression  equal  to  zero.  Solving  for 

cDmin  yiclds * 


W2 


CDmin  ij  -n — jr 

Hrv07*  A 


(4?) 


where  VQ  is  the  speed  for  minimum  sink  rate  (minimum  power). 
The  minimum  power  required  for  level  flight  may  now  be 
written ; 


Po  = 


]/■  + w2 

T~r  VoA  P V0  A 


(43) 
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K-  «-»-.  d»“'d  * tM  *lrCr‘ft  th 
prod»c«  tft.  sin*  «*•  -»  " W>  ^ * 

.ini—  .1-*  »*  ™*  *“  •It”M*d' 


R/S, 


4w 


(44) 


3T  V0A 

Mans  (4?)  and  (44)  into  equation 
Substituting  equations 

.imolified,  more  useful  form  of  the 
(41)  produces  a simple 

power  required  equation; 

. „ 2 (45) 

p,  . W V3  R/so  ♦ 


n vT 


4 V 


ot  eonr...  «-  W <»Ur  «-*  * * f" 

v.i»oity  "*y  * * th*  foU”'1"8  5t"5' 

(46a) 

na9eC^  (46b) 

tan^  * VA/g 

(46c) 

n »\j  (Viv/g)  + 1 

._d  may  be  written  in  terms  of 
Now  the  power  required  may 

the  turning  rate. 

Pr  «=  W R/S0  (V/V0)3  + 

+ 3 VJ  R/S0  Vc 
— TTV  ’ 


The  energy 

fr/ii- 

Er  *{  Pi*  dt 


required  for  the  turn  is* 


Fr  - WTEfR/S0  (V/V0^  + 3 R/S„  (V0/V)) 

f i.  \ 4 4 y 


WY-fL 


3 R/S0  Vn  V 


The  theoretical  turning  rate  for  minimum  energy 
expenditure  may  he  calculated  from  equation  (47)  to 


_L2_ (minimum  energy)  = g|  ,'( V/V0) ^ + 3 (V./V) 


This  theoretical  result,  however,  produces  unaccept- 
ably large  hank  angles  and  cannot  practically  he  used  as 
a turn  criterion.  Tn  the  numerical  analysis  which  follows 
this  section  a hank  angle  of  forty-five  degrees  will  he 
assumed.  Obviously  it  is  easier  for  the  pilot  to  con- 
form to  a constant  hank  angle  than  for  a specified  turn- 
ing rate. 

At  the  high  speed  end  the  pilot  may  wish  to  choose 
any  convenient  speed  less  than  the  aircraft’s  rated  max- 
imum sDeed . A suggested  criteria  for  Vmax  would  be  the 
maneuvering  speed,  or  the  speed  recommended  by  the  man- 
ufacturer for  mild  aerobatics. 

Tn  choosing  the  minimum  velocity  for  the  racetrack 
schedule  the  pilot  may  wish  to  utilize  the  following  tech- 
nique: First,  a desireable  hank  angle  must  be  chosen  and 

equation  (46a)  employed  to  find  the  load  factor  correspond- 
ing to  that  bank  angle.  The  stall  speed  in  a turn  varies 
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directly  with  the  square  root  of  the  load  factor  in  the 


turn? 


V( stall  in  turn)  = V(stall)^"rT 


The  load  factor  in  a 45®  banked  turn  is  n a 1,4. 

For  the  numerical  analysis  Vnin  will  be  taken  to  be  the 
stall  speed  times  the  square  root  of  1 .4.  Since  the  turn 
for  a sailplane  cannot  be  truly  level,  the  load  factor 
will  actually  be  slightly  smaller  and  therefore  a margin 
of  safety  need  not  be  utilised. 

Now  the  power  required  to  complete  the  phugoid 
oscillation  will  be  examined.  From  equation  (43)  with 
the  load  factor  equal  to  one, 

Pr  a R/Sc  W (V3/V3  + 3VC/V)  (50 


The  velocity  for  this  portion  of  the  flight  is  given 
by  equation  (39).  rIpon  substitution  and  expansion,  equa- 
tion (50)  becomes? 

Pr  = R/Sp  W (V/VQ)3  + yPw  eoswt  + 3V&V2  cos^wt 
~~n  L 2vg  4V3 

• + l(4V/VQ)3cos3wt  + 3 V_(V  + ^AVcoswt)  ^ 

■5 


The  energy  required  is  obtained  by  integrating 
equation  (51)  over  one  phugoid  period. 


2lr/w 
I Pr  dt 


Ml  odd  powers  of  cosine  drop  out  in  the  integration. 
Therefore  the  result  of  the  integration  is; 


F>  = R/S0  W1TV  [ V2  + 3(AV/2)2]  + (5 2) 

~ — rr-  ? — 


2 w v. 


2wPn**Vn 


The  total  energy  required  for  the  racetrack  flight 
schedule  is  the  sun  of  equations  (47)  and  (52). 
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The  power  available  from  the  wind  is; 

Pa  - m dVw  V - W dVw  dh  V ( 54) 

dt  dh  dt 


Taking  time  zero  to  be  at  the  beginning  of  the  ascent 
(V(0)  <=  Vmax),  the  equations  for  the  flight  parameters  are 
merely  sinusoidal  functions  of  the  phugoid  frequency. 

V = V + jAVcoswt  (39) 

# = *maxsinwt  (55) 

dh/dt  = (dh/dt)maxsin>rt  (5^) 

Evaluation  of  (dh/dt)  . is  similar  as  for  the  cir- 

ma* 

cular  case. 


T'/tt J 

Ah  = J (dh/dt)maxsinwt  dt 
•/o 

w 

but  mg  Ah  = ^(Vm^x  - Vm^n  ) 

so  (dh/dt)way  = - W(Vm;x  - Vmfn  ) (57) 


Substituting  equations  (39),  (5^),  and  (57)  into 
equation  (54)  and  assuming  a wind  shear  of  constant 
value  k in  the  vicinity  of  the  flight  path; 

Pa  - -Wkw(Vnax  - Vmfn)  sinwt  (V  + f AVcoswt) 


(5B) 


The  negative  sign  indicates  that  power  is  being  sup- 
plied to  the  aircraft  by  the  wind.  The  energy  provided  by 
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the  climb  is  the  same  as  that  provided  by  the  glide. 


Therefore  the  total  energy  supplied  by  the  wind  along  the 


straight  portions  of  the  flight  schedule  isi 


Ea  = 2(  - Pa  dt 


Ea  - 4 W k V2 
F2 


In  the  level  turns  the  wind  will  supply  or  extract 


energy  from  the  aircraft.  The  power  contribution  of  the 


wind  is  found  from; 


Pa  = W dVg  y 


8 dt“ 

In  the  level  turn 


Vg  = V + VsosV 


V = constant 


-Vw  sin'*'  dY_ 


so  Pa  * - V Vtf  WjvsinV 


For  the  turn  at  the  high  speed  end  (turning  into  the 
wind)  the  energy  provided  by  the  wind  ist 


-Jpa  dt  = - 1 jpa  dV 


- _ 2 V V W 
max  vw  w 


Equation  (6l)  shows  that  energy  is  extracted  from  the 


aircraft  in  the  high  speed  turn  into  the  wind.  For  the 
turn  at  low  speed  away  from  the  wind  the  integration 
shows  that  energy  is  supplied  to  the  aircraft. 


^ vw  w 


(62) 


Combining  equations  (6l)  and  (62)  yields  the  total 
energy  contribution  in  the  turns  due  to  the  wind; 

Ea  = ? W AV  Vw  (fi3) 

Thus  the  total  energy  available  from  the  racetrack 
flight  path  is  found  by  adding  equations  (59)  and  (6?); 


Ea  = 4 W A V / V2. 

-2T-ik  + v- 


(64) 


' / 

Comparing  equation  (64)  with  equation  (38)  shows 
that  the  racetrack  flight  schedule  provides  more  energy 
from  the  wind  by  a factor  of  4/ IT  . 


EXPERIMENTAL  RESULTS  AND  NUMERICAL  ANALYSIS 


Experimental  investigation  was  conducted  in  a Beech- 
craft  Bonanza  Model  F-33A  airplane.  This  is  a four  seat 
low  wing  general  aviation  aircraft.  It  is  powered  by  a 
Continental  I0-520-BA  reciprocating  engine  with  rated  max- 
imum horsepower  of  285  «t  sea  level.  The  propeller  is  a 
constant  speed  two  bladed  type.  All  flights  for  data  were 
made  with  gear  and  flaps  retracted. 

The  aircraft  was  equipped  with  a drift  sight  located 
between  the  front  and  real*  seats  on  the  pilot's  side.  The 
drift  sight  was  used  to  measure  drift  angles  in  straight 
and  level  flight.  By  triangulation  methods  wind  speedsiat 
flight  altitudes  were  determined. 

Figure  7 shows  an  upper  level  shear  discovered  on 
24  March  1976.  A non-linear  shear  interpretation  has  been 
indicated.  The  basis  for  this  representation  of  the  bhear 
profile  is  the  results  indicated  in  Figure  3 of  the  Intro- 
duction, The  presence  of  such  a non-linear  shear  is  also 
indicated  by  Figure  8. 

Figure  8 depicts  experimental  results  from  glides  in 
the  area  of  the  24  March  shear.  The  glides  were  made  power 
off  both  with  the  wind  and  against  the  wind  to  verify  the 
effect  of  wind  shear  predicted  by  equation  (17).  Altitude 
and  elapsed  time  were  recorded  continuously  throughout  the 


iirure  7.  Upner  lev«  1 shear  measured  2 k March  1 a7^. 
Linear  shear  value  .003^9. 


FIGURE  8.  SINK  RATE  VERSUS  SHEAR 


pi ides.  Sink  rates  at  various  altitudes  (i.e.  shear  values) 


have  been  plotted  in  linear  fern  against  the  shears  oorres 


ponding  to  those  altitudes 


The  plot  for  the  glide  with  the  wind  indicates  that  sink 


rate  is  reduced  with  an  increased  shear.  Also  quite  evident 


in  Figure  8 is  difference  between  gliding  with  the  wind  and 


On  14  January  10?8  drift  sight  readings  were  taken  to 


verify  the  existence  of  a boundary  layer  shear.  Table  ?. 


shows  the  wind  velocities  measured  over  the  Chesapeake  Bay 


Boundary  Layer  Wind  Velocities  Measured  14  Jan.  197^ 


Altitude  feet  Wind  Velocity  ft/sec 


The  strength  coefficient  for  Davenport's  power  equation 


was  evaluated  at  each  altitude  and  averaged  to  discover  the 


equation  for  the  boundary  layer 
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Figure  Q is  a plot  of  equation  (65a)  with  data  points 
indicated.  Correlation  is  good  and  Davenport's  power 
equation  seems  to  adequately  describe  the  boundary  layer 
shear  profile.  Table  3 lists  boundary  layer  wind  parameters 
taken  from  equations  (65a)  and  (65b), 

Table  3 

Wind  Parameters  In  The  Boundary  Layer 
Altitude  feet  Wind  Velocity  ft/sec  Shear  sec-1  ! 


100 

32.2 

.046 

200 

35.6 

.025 

300 

37.7 

.018 

400 

39.3 

.014 

500 

40.5 

.012 

Utilizing  the  boundary  layer  data  obtained  on  14  Jan- 
uary for  wind  parameters,  the  Beechcraft  Bonanza  was  ev- 
aluated for  dynamic  soaring  performance  using  the  race- 
track flight  technique, 

A computer  program  was  used  to  solve  equations  (53) 
and  (64)  for  the  aircraft  flying  in  the  region  described 
by  Table  3 (i.e,  shears  from  approximately  ,01  to  .08), 

In  addition,  data  from  Appendix  1 was  used;  the  air- 
plane's performance  at  simulated  superior  glide  ratios  was 
investigated.  Tables  4 and  5 indicate  parameters  used  for 
the  numerical  analysis. 
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1 Table  4 

i Beeehcraft  Bonanza  F-33A  Performance  Data 

! Stall  Speed,  Power  Off 

74  MPH 

Design  Maneuvering  Speed 

152  MPH 

Phugoid  Period 

23.7  seconds 

Vmax 

152  MPH 

^min 

88  MPH 

V 

120  MPH 

A V 

32  MPH 

Table  5 

Beeehcraft  Bonanza  F-33A  Sink  Rate  Data 
At  Various  Power  Settings 
(From  Appendix  1) 

Power  Setting 

R/S„  ft/s 

MPH 

0* 

14.3 

88 

15* 

l 

to 

* 

20* 

i 

8.8 

)0 

OO 

25* 

6.2 

98 

30* 

4.6 

99 

35* 

3.0 

99 

40* 

1.5 

101 

Figure  10  is  a plot  of  altitude  gained  (energy  avail- 


able divided  by  weight)  in  one  racetrack  cycle.  The  alti- 
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tude  gained  due  to  wind  shear  is  independant  of  the  sink 
rate  performance  of  the  aircraft  hut  rather  solely  depends 
upon  wind  characteristics  and  velocity  parameters  chosen 


! 


i 

i 


i 

- 

i 


by  the  pilot. 

Figure  l1  shows  the  result  of  subtracting  the  alti- 
tude lost  from  the  altitude  gained  for  one  racetrack  cycle. 
As  is  clearly  evident,  at  low  power  settings  (low  glide 
ratios)  the  net  loss  in  altitude  is  totally  unacceptable. 
Only  with  extremely  high  simulated  glide  ratios  does  the 
altitude  lost  begin  to  diminish  to  an  acceptable  value. 
Sailplanes  with  glide  ratios  of  the  order  of  magnitude 
necessary  to  produce  a net  increase  in  altitude  per  cycle 
do  not  exist  today.  Further  numerical  analysis  of  the 
dynamic  soaring  problem,  utilizing  actual  sailplane  param- 
eters is  warranted,  however,  in  order  to  determine  the 
true  effect  of  the  racetrack  schedule  upon  high  performance 
sailplanes. 
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CONCLUSIONS 


1.  Wind  shear  nay  be  used  to  increase  the  perform- 
ance of  an  aircraft.  This  increase  in  performance  may 
be  expressed  mathematically,  and  experimental  observations 
concur  with  the  theoretical  predictions, 

2,  The  wind  shear  found  in  the  upper  atmosphere 
(upper  level  shear)  is  not  strong  enough  to  provide  the 
energy  required  to  maintain  continuous  non-powered  flight 
for  even  the  best  high  performance  aircraft.  This  does 
not  rule  out  the  possibility  of  a jet  stream  shear,  how- 
ever, which  was  not  examined  in  this  report. 

3*  Wind  shear  found  in  the  lower  part  of  the  boundary 
layer  is  of  the  order  of  magnitude  required  to  maintain  con- 
tinuous non-powered  flight, 

A racetrack  flight  pattern  was  found  to  be  more 
condusive  to  practical  application  than  the  circular  case 
from  two  standpoints.  First,  it  turned  out  to  be  more 
efficient  theoretically.  Secondly  it  was  much  easier  for 
the  pilot  to  conform  to  the  theoretical  model  for  the  race- 
track than  the  circle, 

5.  The  successful  employment  of  the  racetrack  flight 
technique  described  in  this  paper  greatly  depends  upon  the 
aircraft’s  performance  parameters.  High  glide  ratio  sail- 
planes would  definitely  be  a requirement  for  achieving  the 
objective  of  this  study. 
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The  numerical  analysis  of  the  racetrack  flight 
schedule  suggests  a possible  method  for  increasing  the  on 
station  time  of  patrol  aircraft  operating  in  wind  shear 
weather  conditions.  Py  properly  employing  wind  shear  at 
a reduced  power  setting  a reduction  in  fuel  consumption 
may  be  achieved. 

7.  The  possibility  of  utilizing  vertical  wind  currents 
(updrafts)  in  conjunction  with  horizontal  wind  shear  was 
not  explored  by  this  paper.  Such  a topic  is  suggested  as 
a basis  for  further  investigation. 
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APPENDIX  1 


I 

I \ 


SIMULATION  OF  SUPERIOR  GLIDE  RATIOS 

I 

Flight  testing  the  concept  of  dynamic  soaring  in  wind 
shear  with  a sailplane  would  require  a great  deal  of  time 
and  expense.  In  addition,  utilizing  one  sailplane  limits 
the  investigator  to  a single  set  of  aircraft  performance 
parameters.  It  is  desirable  to  employ  a single  aircraft 
. with  variable  sailplane  characteristics  which  can  make 
repeated  test  runs  in  a single  flight.  It  is  possible  to 
do  Just  that  with  a powered  aircraft.  Lift  to  drag  ratios 
exceeding  present  sailplane  maximum*  may  be  reached  by 
applying  a fraction  of  the  power  necessary  to  maintain 
level  flight. 

The  aircraft  to  be  used  for,  experimenraj.  tesxing  of 
dynamic  soaring  analytic  results  was  the  United  States 
Naval  Academy's  Eeechcraft  Bonanza.  A brief  description 
of  the  process  used  to  calibrate  various  aircraft  power 
settings  with  glide  ratios  and  results  are  contained  herein. 
Actual  intermediate  curves  and  computer  printouts  are 
available  from  the  author  in  a seperate  unpublished  report 
sr.titled  "Simulation  of  Superior  Lift  To  Drag  Character- 
istics With  A Powered  Aircraft", 

Power  required  and  power  available  curves  for  the 
Beechcraft  Bonanza  were  determined  from  actual  inflight 
testing  on  October  9,  1975*  Power  required  was  ascertained 
from  a linear  plot  of  V**  versus  P*V.  Power  available  was 
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found  from  an  acceleration  run  which  produced  excess 
power  data.  Excess  power,  when  added  to  power  required, 
produces  the  power  available  curve. 

Percentage  curves  of  the  power  available  curve  were 
then  computed.  These  curves  then  became  the  zero  line 
for  power  required  and  points  from  the  former  were  sub- 
tracted from  the  latter.  The  result  was  a reduction  in 
the  effective  power  required  and  a resultant  increase  in 
glide  ratio.  Figure  6 displays  effective  power  required 
curves  (multiplied  by  550  horse power /ft- lb  and  divided  by 
standard  Beechcraft  weight  of  3400  pounds  to  yield  sink 
rates)  at  power  settings  of  various  levels.  Results  of 
uaictu  a i. j ons  are  Tabulated  lr.  Table  1. 


TABLE  1. 
Glide  Ratios  At  Various 

Power  Settings 

Power  Setting 

L/Dmax 

o$ 

11.5? 

15* 

15.28 

20$ 

18.89 

25$ 

23.66 

30$ 

31.27 

35$ 

48.37 

40$ 

• 1 

96.73 
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APPENDIX  2 
NOMENCLATURE 

a exponent  of  boundary  layer  shear  equation?  based 
upon  terrain  considerations. 

ax  acceleration  in  the  x direction. 

A ft*  e 

* aspect  ratio. 

cDmin  dr*E  coefficient  at  eero  lift. 

CTi  lift  coefficient, 

D drag. 

e span  efficiency  factor. 

Ea  energy  available. 

Er  energy  required. 

E descent  efficiency  factor  due  to  wind  shear. 

F force  acting  on  the  aircraft  due  to  wind  shear. 

Fx  force  in  the  x direction. 

Fy  force  in  the  y direction, 

f phugoid  cyclic  frequency, 

P acceleration  due  to  gravity, 

h altitude. 

k strength  coefficient  of  boundary  layer  shear  power 

equation. 

KE  kinetic  energy. 

L lift. 

m,M  mass, 

n load  factor. 
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P power . 

Pa  power  available. 

Pr  power  required, 

q dynamic  pressure. 

R/D  rate  of  descent. 

R/Sq  minimum  rate  of  sink. 

S wing  surface  area, 

t time, 

V true  airspeed. 

Vg  ground  speed. 

V0  speed  for  minimum  sink  rate  (minimum  power), 
Vw  wind  speed, 

Vy  vertical  velocity, 

w phugoid  angular  frequency. 

W weight. 

i altitude, 

flight  angle , 

r ips 

p density. 

flight  heading  with  respect  to  wind  direction. 
XL  angular  rotation  rate  (turning  rate). 

<6  bank  angle. 
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